To appear in the Astrophysical Journal 



Magnetic Field Structure of the HH 1-2 Region: Near-Infrared Polarimetry of Point-Like Sources 
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ABSTRACT 

The HH 1-2 region in the L1641 molecular cloud was observed in the near-IR J, H, and K s bands, and imaging 
polarimetry was performed. Seventy six point-like sources were detected in all three bands. The near-IR polarizations of 
these sources seem to be caused mostly by the dichroic extinction. Using a color-color diagram, reddened sources with 
little infrared excess were selected to trace the magnetic field structure of the molecular cloud. The mean polarization 
position angle of these sources is about 111 , which is interpreted as the projected direction of the magnetic field in 
the observed region of the cloud. The distribution of the polarization angle has a dispersion of about 11°, which is 
smaller than what was measured in previous studies. This small dispersion gives a rough estimate of the strength of 
the magnetic field to be about 130 /iG and suggests that the global magnetic field in this region is quite regular and 
straight. In contrast, the outflows driven by young stellar objects in this region seem to have no preferred orientation. 
This discrepancy suggests that the magnetic field in the LI 641 molecular cloud does not dictate the orientation of the 
protostars forming inside. 

Subject headings: infrared: stars — ISM: individual (HH 1-2) — ISM: structure — polarization — stars: formation 



1. INTRODUCTION 

Magnetic fields play a crucial role in various astrophysi- 
cal processes, including the evolution of interstellar molecular 
clouds and star formation (Shu et al. 1987; Bergin & Tafalla 
2007; McKee & Ostriker 2007). One of the problems related 
to star formation concerns the competition between magnetic 
and turbulent forces (Mac Low & Klessen 2004). The magnetic 
field direction can be measured by observing the dichroic po- 
larization of background stars in the optical and near-IR bands 
and/or the linearly polarized emission from the dust grains in the 
mid-IR and far-IR bands (Davis & Greenstein 1951; Matthews 
& Wilson 2000). The large-scale alignment of dust grains with 
the magnetic field is known to be the cause of the dichroic ex- 
tinction and the interstellar polarization seen in the direction 
of background sources. Because of the low extinction, near-IR 
imaging polarimetry is particularly useful in tracing the dichroic 
polarization of background stars and embedded sources seen 
through dense clouds (Vrba et al. 1976; Wilking et al. 1979; 
Tamura et al. 1987; Kandori et al. 2007). Since both dichroic 
extinction and scattering processes can contribute to the polar- 
ization of embedded sources, multiwavelength polarimetry can 
be useful in discriminating between the two mechanisms (Casali 
1995). 

The LI 641 cloud is one of the nearest giant molecular clouds 
and is a site of active star formation (Kutner et al. 1977; Mad- 
dalena et al. 1986; Strom et al. 1989; Morgan & Bally 1991; 
Sakamoto et al. 1997; Zavagno et al. 1997). The role of mag- 
netic field in the star formation activity of LI 641 is complicated. 
With visual polarimetry of background stars, Vrba et al. (1988) 
found that the dispersion in position angles is large (33°) and 
suggested that the role of magnetic field in the global scale is 
only incidental. However, they also found that the outflows in 
L1641 tend to be parallel to the field direction and suggested 
that the role of magnetic field is important in the local scale. In 
contrast, Casali (1995) performed near-IR polarimetry of young 



stellar objects (YSOs) in LI 641 and found that the alignment of 
polarization vectors is poor, which suggests that the magnetic 
field was not dominant in the collapse dynamics. To understand 
the situation better, it was suggested that more extensive po- 
larimetry in the vicinity of each outflow and YSO is necessary 
(Vrba etal. 1988). 

One of the well studied parts of the LI 641 cloud is the re- 
gion around the reflection nebula NGC 1999 and the Herbig- 
Haro objects HH 1-2. This region contains several YSOs and 
outflows (Herbig 1951; Haro 1952; Warren-Smith et al. 1980; 
Strom et al. 1989; Corcoran & Ray 1995; Choi & Zhou 1997; 
Rodriguez et al. 2000). The magnetic field structure in the 
HH 1-2 region has been studied based on optical polarizations 
of point sources (Strom et al. 1985; Warren-Smith & Scarrott 
1999, hereafter WS). They found that the local magnetic field is 
directed roughly along the axis of the HH 1-2 outflow, but the 
number of detectable stars was too small because of the large 
obscuration in this region. 

In this paper, we present a wide-field near-IR polarimetry of 
the HH 1-2 region. In Section 2 we describe the observations 
and data reduction. In Section 3 we present the results of the 
polarimetry of point-like sources. In Section 4 we discuss the 
magnetic field structure and the star-forming activity in the HH 
1-2 region. A summary is given in Section 5. 

2. OBSERVATIONS 

The observations toward the HH 1-2 region were carried out 
using the SIRPOL imaging polarimeter on the Infrared Survey 
Facility (IRSF) 1.4 m telescope at the South African Astronom- 
ical Observatory. SIRPOL consists of a single-beam polarime- 
ter (an achromatic half-wave plate rotator unit and a polarizer) 
and an imaging camera (Nagayama et al. 2003). The camera, 
SIRIUS, has three 1024 x 1024 HgCdTe infrared detectors. 
IRSF/SIRPOL enables deep and wide-field (7.'7 x T.7 with a 
scale of 0/!45 pixeU 1 ) imaging polarimetry at the J, H, and K s 
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[See http://minho.kasi.re.kr/Publications.html for the original high-quality figure.] 



Fig. 1. — Color composite Stokes / image of the HH 1-2 region in the J(blue), H(green), and K s (red) bands from the IRSF/SIRPOL 
observations. 



bands simultaneously (Kandori et al. 2006). 

The observations were made on the night of 2008 January 
9. We performed 20 s exposures at 4 wave-plate angles (in the 
sequence of 0°, 45°, 22? 5, and 67? 5) at 10 dithered positions 
for each set. The same observation sets were repeated 10 times 
toward the target object and the sky backgrounds for a better 
signal-to-noise ratio. The total integration time was 2000 s per 
wave plate angle. The typical seeing size during the observa- 
tions was ~1"3 in the J band. The polarization efficiencies of 
SIRPOL are stable over several years, and the instrumental po- 
larization is negligible (Kandori et al. 2006). The efficiencies 
were measured in 2007 December during a maintenance period, 
just a few days before our observing run, and were the same as 



the values reported by Kandori et al. (2006). 

The data were processed using IRAF in the same manner 
as described by Kandori et al. (2006), which included dark- 
field subtraction, flat-field correction, median sky subtraction, 
and frame registration. Figure 1 shows the J-H-K s color com- 
posite intensity image of the 8' x 8' region around HH 1-2 
(hereafter the HH 1-2 field). Many point-like sources in the HH 
1-2 region were detected. In addition to the point-like sources, 
we detected HH 1-2, NGC 1999, and other nebulosities, but the 
polarimetric studies of these extended sources will be presented 
elsewhere in the future. 




3. RESULTS 

3.1. Photometry 

The IRAF DAOPHOT package was used for source detec- 
tion and photometry (Stetson 1987). The DAOPHOT program 
automatically detected point-like sources with peak intensities 
greater than 10a above the local sky background, where a is the 
rms uncertainty. The automatic detection procedure misidenti- 
fied some spurious sources and missed some real sources, and 
the source list was corrected by visually inspecting the im- 
ages carefully. The pixel coordinates of the detected sources 
were matched with the celestial coordinates of their counter- 
parts in the Two Micron All Sky Survey (2MASS) Point Source 



Catalog. The IRAF IMCOORDS package was applied to the 
matched list to obtain plate transform parameters. The rms un- 
certainty in the coordinate transformation was ~0l! 1 . 

Aperture photometry was performed again with the resulting 
images. The aperture radius was 3 pixels, and the sky annulus 
was set to 10 pixels with a 5 pixel width. The resulting list con- 
tains 76 sources whose photometric uncertainties are less than 
0.1 mag in all three bands (Table 1). These point-like sources 
are labeled in Figure 2. Four bright sources (V380 Ori, the C-S 
star, N 3 SK 50, and an unnamed star ~10" south of source 26) 
were saturated, and they were excluded from the list. 

The Stokes I intensity of each point-like source was calcu- 
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lated by 



I — 2 + ^22.5 + lib + ^67.5) 



(1) 



where I a is the intensity with the half wave plate oriented at a . 
The magnitude and color of the photometry were transformed 
into the 2MASS system by 

MAG 2M ass = MAGirsf + <*i x COLORi RSF + ft (2) 



and 



COLOR2MASS = ol 2 x COLORirsf + fa, 



(3) 



where MAGirsf is the instrumental magnitude from the IRSF 
images, and MAG2MASS is the magnitude from the 2MASS 
Point Source Catalog. The parameters were determined by fit- 
ting the data using a robust least absolute deviation method. 
For the magnitudes, a x = 0.017, -0.064, and 0.001, and ft = 
--4.986, -4.717, and -5.375 for J, H, and K s , respectively. For 
the colors, a 2 = 1.007 and 0.960, and fa = -0.261 and 0.664 for 
J — H and H — K s , respectively. The coefficients ft and fa 
include both the zero point correction and aperture correction. 
The derived magnitudes are listed in Table 1. The 10ct limiting 
magnitudes were 19.6, 18.7, and 17.3 for J, H, and K s , respec- 
tively. 

3.2. Polarimetry 

Aperture polarimetry was carried out on the combined in- 
tensity images for each wave plate angle, instead of using the 
Stokes Q and U images. This is because the center of the 
sources cannot be determined satisfactorily on the Q and U im- 
ages. From the aperture photometries on each wave plate angle 
image, the Stokes parameters of each point-like source were de- 
rived by 

Q = Io- h 5 (4) 



and 



U = l 



22.5 



'67. 5- 



(5) 



The aperture and sky radius were the same as those used in the 
photometry of I images. The degree of polarization, P, and the 
polarization position angle, 9, can be calculated by 



Pn = 



and 



VQ 2 + u 2 



P 2 - 5P\ 



1 u 

-arctan-, 



(6) 



(7) 



(8) 



where 5P is the uncertainty in P . Equation (7) is necessary 
to debias the polarization degree (Wardle & Kronberg 1974). 
Finally, P was corrected using the polarization efficiencies of 
SIRPOL: 95.5 %, 96.3 %, and 98.5 % at J, H, and K s , respec- 
tively (Kandori et al. 2006). 

Table 2 shows the derived source parameters. The uncertain- 
ties given in Table 2 (and elsewhere in this paper) are la val- 
ues. Figures 3-5 show the polarization vector maps of point-like 
sources superposed on the I images. For the sources with PI8P 
> 4 and P < 9% (21 sources), the correlation coefficients are 
0.91 for (9 H , 9k s ) and 0.97 for (6 H , 6j). Figure 6 shows the his- 
tograms and Gaussian fits for the polarization position angles. 
Each of the three histograms shows a single peak at ~111°. 
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[See http://minho.kasi.re.kr/Publications.html for the original high-quality figure.] 

Fig. 3. — Stokes / image of the J band with polarization vec- 
tors. The length of the vectors is proportional to the degree of 
polarization. Shown in the upper right corner is a 10% vector. 
Note that there are bad pixel clusters around the upper-left and 
upper-right corners and the middle of the right boundary. 
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Fig. 4. — The same as Fig. 3 for the H band. 



Note that the dispersion in 9 is smallest in the H band. In addi- 
tion, for most sources, the signal-to-noise ratio (P/5P) is higher 
in the H band than the other bands. The H band polarimetry is 
more reliable than those of the other bands probably because 
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[See http://minho.kasi.re.kr/Publications.html for the original high-quality figure.] 

Fig. 5. — The same as Fig. 3 for the K s band. 
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Fig. 6. — Histograms of polarization position angles for the J, 
H, and K s bands. All the sources in Table 2 with detected po- 
larization are included. Dot-dashed curves: Gaussian fits. The 
peak angle and dispersion are 1 12° and 16° for J, 111 and 13° 
for H, and 1 10° and 19° for K s . 

the contamination from extended nebulosity is smaller in the H 
band than in the J band and because the dichroic polarization 
is more efficient in the H band than in the K s band. Therefore, 
our discussion in Section 4 will be mainly based on the H band 
data. 

The relation between the polarimetric and spectral data may 
be useful in understanding the nature of polarization. The de- 
gree of polarization appears to be correlated with near-IR colors 
(Fig. 7). The empirical relation for the upper limit of interstellar 
polarization suggested by Jones (1989) is 

i^max =tanh jl.5£(if-ir)^jj, (9) 

where r\ = 0.875 and E(H — K) is the reddening owing to ex- 
tinction. Most of the sources are within this limit (Fig. 7). A 
few sources are above the P max limit, but their uncertainties are 
large. The near-IR polarization-to-extinction efficiency of the 
point-like sources in the HH 1-2 field is consistent with that 
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Fig. 7. — Degree of polarization vs. H — K s color. Dot-dashed 
lines: empirical upper limits (P max ; Jones 1989). 

caused by aligned dust grains in the dense interstellar medium. 
Therefore, their polarizations are likely dominated by the inter- 
stellar dichroic extinction, and the intrinsic polarization, if any, 
did not significantly enhance the degree of polarization. How- 
ever, this result does not completely exclude the possibility that 
some of the sources have intrinsic polarization because depolar- 
ization is also possible. 

4. DISCUSSION 

4.1. Comparison with Previous Studies 

Polarimetry of bright point-like sources in the HH 1-2 re- 
gion was reported previously by several authors. These studies 
covered larger regions than our study, but they were much shal- 
lower. Strom et al. (1985) carried out I band polarimetry and 
measured the polarization position angle of two sources: 157° 
for the C-S star and 135° for N 3 SK 50. Casali (1995) measured 
the K band polarization angle of two sources: 105° for V380 
Ori and 154° for N 3 SK 50. These three sources were saturated 
in our observations, and no direct comparison is possible. 

WS presented a broad-band (450-1000 nm) polarimetry of 
eight point-like sources in a larger (~10') region. The polariza- 
tion position angle averaged about 130° with a dispersion of 
30°. The polarization angle of bright sources were 125° for the 
C-S star and 128° for N 3 SK 50. Direct comparisons for polar- 
ization angles are possible for three sources, while meaningful 
comparisons for polarization degrees are difficult due to wave- 
length dependence. WS 1 (source 1 of WS) corresponds to our 
source 8, and the polarization position angle of WS (145° ±5°) 
is somewhat larger than our near-IR measurements (95-121°). 
WS 5 and WS 6 correspond to our sources 10 and 58, respec- 
tively, and their polarization orientations agree reasonably well, 
within 5°. It is not clear what caused the difference of source 
8. One of the possibilities is that the polarization of source 8 
is caused mostly by the scattering process rather than dichroic 
extinction, as this source shows little extinction (see Section 4.7 
for more discussions). 

Based on the broad-band polarimetry, WS suggested that the 
magnetic field in the HH 1-2 region is oriented at a position 
angle of 126°. This interpretation, however, should be corrob- 
orated by deeper observations because the size of their source 
sample was too small. Our observations can provide statistically 
more significant interpretations, which are presented below. 

4.2. Source Classification 

To study the magnetic field structure of molecular clouds 
through the interstellar polarization caused by dichroic extinc- 
tion, it is necessary to select sources without intrinsic polar- 
ization. YSOs in the cloud can exhibit a substantial degree of 
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intrinsic polarization caused by circumstellar material. Such 
sources may show a large amount of infrared excess emission. 
Therefore, it is important to classify sources, and the multiwave- 
length photometry can be useful. 

Figure 8 shows a color-color diagram for all the sources de- 
tected in all three bands. The diagram was divided into several 
domains. Based on the location in this diagram, sources can be 
classified into a few groups (Lada & Adams 1992). The area 
near the locus of main-sequence/giant stars is called domain 
AO. Sources in domain AO are either field stars (dwarfs and gi- 
ants) or pre-main- sequence (PMS) stars with little infrared ex- 
cess (weak-lined T Tauri stars and some classical T Tauri stars) 
and with little reddening. There is a clear gap just above domain 
AO, and the area above this gap in the direction of the redden- 
ing vector is called domain Ar. Sources in domain Ar are either 
field stars or PMS stars with little infrared excess and with sub- 
stantial reddening. Domain B is the area next to domain Ar in 
the direction of higher H — K s (to the right) and above the lo- 
cus of classical T Tauri stars. Sources in domain B are PMS 
stars with infrared excess emission from disks. Domain C is the 
area next to domain B to the right. Sources in domain C are 
infrared protostars or Class I sources. Herbig AeBe stars tend 
to occupy lower parts of domains B and C. This classification 
based on the color-color diagram, however, is far from perfect. 
A certain fraction of classical T Tauri stars may reside in do- 
mains AO or Ar, some protostars can be found in domain B, and 
some extremely reddened AeBe stars may be found among pro- 
tostars (Lada & Adams 1992). This "contamination" will even- 
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Fig. 8. — Color-color diagram of the point-like sources in the 
HH 1-2 field. Filled circles: point-like sources from this work 
(Table 1). Open circles: bright sources (Strom et al. 1989; Car- 
penter et al. 2001). Solid curve: locus of main sequence and 
giant branch stars (Bessell & Brett 1988). Dotted line: locus of 
classical T Tauri stars (Meyer et al. 1997). Dashed lines: bound- 
aries between domains Ar, B, and C (see Section 4.2). Solid 
line: reddening vector. 



tually contribute to the uncertainty in statistical quantities de- 
rived from the classification, but the estimation of this uncer- 
tainty is beyond the scope of this paper. 

There are thirteen sources in domain AO, and they are col- 
lectively called group AO. They are either foreground stars or 
those seen along lines of sight with little extinction. 

Fifty seven sources were found in domain Ar. These sources 
(group Ar) are either background stars or PMS stars in the 
LI 641 cloud. They are the most useful sources for the study 
of the magnetic fields in the cloud (Section 4.4). 

Six sources were found in domain B. These sources (group 
B) may be PMS stars associated with the L1641 cloud. Source 
10 (WS 5) is the emission-line star AY Ori (Wouterloot & Brand 
1992) and also source 13932 of Carpenter et al. (2001). Source 
21 is source 14843 of Carpenter et al. (2001). In addition, two 
of the brightest objects in this field, the C-S star and N 3 SK 50, 
also belong to group B. 

None of our sources are located in domain C. However, this 
nondetection does not mean that there is no protostar in the 
HH 1-2 field. Some protostars (for example, HH 1-2 VLA 1) 
are deeply embedded and undetectable in the near-IR bands. In 
addition, the red protostars tend to be missed at shorter wave- 
lengths (J and/or H bands). The Herbig Ae star V380 Ori is 
located in the lower part of domain C, as expected. 

4.3. Source of Polarization 

To study the magnetic field structure of LI 641, we will use 
the distribution of the polarization of group Ar sources (Sec- 
tion 4.4). To do that, the source of polarization needs some dis- 
cussion. The distribution of polarization angle (Fig. 6) shows a 
well-defined single peak, which suggests that the cause of polar- 
ization is relatively simple. Especially for the group Ar source, 
the polarization seems to be mostly caused by the LI 641 cloud, 
based on several lines of evidence described below. 

First, optically thick molecular lines (such as CO and H2CO) 
observed toward the HH 1-2 region show a single velocity com- 
ponent at ~8 km s _1 (Loren et al. 1979; Snell & Edwards 
1982), which suggests that LI 641 is the only molecular cloud 
that is optically thick enough to cause the dichroic extinction. 
The spectrum of 13 CO, however, shows that there are two ve- 
locity components: a stronger one at ^8.5 km s _1 and a weaker 
one at ^7 km s _1 (Edwards & Snell 1984), and both compo- 
nents belong to the L1641 cloud (Sakamoto et al. 1997). The 
stronger component is directly related to the dense gas in the 
HH 1-2 region, and the weaker component is probably related 
to the dense gas with a column density peak located ^20' north 
of HH 1-2 (Takaba et al. 1986). Since the distribution of the po- 
larization angle shows a well-defined single peak, one of them 
(most likely the ~8.5 km s^ 1 component) probably dominates 
in the polarization process. Alternatively, the magnetic field di- 
rection of the two components may be very similar. 

Second, the peak polarization angle is insensitive to the 
amount of extinction. If group Ar sources would have a pre- 
ferred direction of polarization even before their IR photons are 
affected by the LI 641 cloud, sources with a small extinction 
would show the effect of this "background" polarization direc- 
tion, while sources with a large extinction would reflect only the 
polarization caused by LI 641 . Figure 9 shows the histograms of 
the polarization position angle for three subgroups of sources 
grouped by the amount of extinction. All the three subgroups 
show a peak angle at ~1 10°. Therefore, the polarization of the 
group Ar sources are mostly caused by the L1641 cloud only, 
and the distribution of the polarization angle may be a very good 
tracer of the magnetic field structure of LI 641. 
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Fig. 9. — Histograms of the H band polarization position an- 
gles of group Ar. The sources were grouped by the amount of 
extinction: Ay = 2-7, 7-11, and 1 1-18. The subgroups contain 
15, 15, and 14 sources, respectively. 



Third, the Galactic latitude of L1641 is high (b = -19?8). It 
is unlikely that light from distant luminous giants would suffer 
a significant amount of dichroic extinction by any other cloud 
behind L1641. 

For some clouds, the polarization direction of background 
stars can be caused by several sources of polarization along the 
line of sight. For example, in the direction of the Southern Coal- 
sack dark cloud (b w -1°), there are at least three components 
of polarization (Andersson & Potter 2005). In such cases, the 
interpretation of the polarization and its relation with the mag- 
netic field structure can be quite complicated. In the case of the 
HH 1-2 region, the polarization seems to be mostly caused by 
LI 641, and the relation between the polarization direction and 
the magnetic field structure is relatively straightforward, as dis- 
cussed in the next section. 

4.4. Magnetic Field Structure 

The sources in group Ar are best for studying the magnetic 
field structure of the molecular cloud because they are subject to 
dichroic extinction and because they would have relatively little 
intrinsic polarization. Figure 10 shows the histogram and Gaus- 
sian fit for the polarization position angles of group Ar sources. 
The peak angle is 111 , and the dispersion is 11°. Comparing 
Figures 6 and 10, it is clear that selecting only group Ar sources 
makes the statistical noise smaller: outliers disappeared, and 
the dispersion became smaller. Therefore, we suggest that the 
global magnetic field in the HH 1-2 region is oriented at a po- 
sition angle of ~1 1 1°. This orientation is consistent with the 
large-scale field structure of L1641 (Vrba et al. 1988). Previous 
studies of the HH 1-2 region (Strom et al. 1985; WS) suggested 
larger position angles and larger dispersions because their sam- 
ple sizes were too small and because they included bright PMS 
stars in the sample. 

To see whether there is a systematic gradient of magnetic 
field orientation over the imaged field, we divided the HH 1- 
2 field into 16 (4 x 4) subregions. In each subregion, group 
Ar sources were selected, and their polarization position angles 
were averaged. Inspection of the resulting distribution of po- 
larization angle (Fig. 11) did not reveal any systematic trend. 
Therefore, we suggest that the measured dispersion of 1 1 rep- 
resents the local variation of magnetic field orientation. Here, 
"local" means each line of sight, though it should be noted that 
the variation along the line of sight would affect the polarization 
in an integrative way. 

Although the polarization measurement does not provide a 
direct estimate of the magnetic field strength at each data point 
in the image, a rough estimation over a large region is pos- 
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Fig. 10. — Histogram of polarization position angles for the 
group Ar sources in the H band. Dot-dashed curve: Gaussian 
fit. 
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Fig. 11. — Average polarization position angle of group Ar 
sources in 2' x 2' subregions of the HH 1-2 field. The num- 
ber of sources in each subregion is labeled. 



sible by statistically comparing the dispersion of polarization 
orientation with the degree of turbulence in the cloud (Chan- 
drasekhar & Fermi 1953). Assuming that velocity perturbations 
are isotropic, the strength of the magnetic field projected on the 
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plane of the sky can be calculated by 



Br, 



(10) 



where Q is a factor to account for various averaging effects, p is 
the mean density of the cloud, <5vi os is the rms line-of-sight ve- 
locity, and 89 is the dispersion of polarization angles. Ostriker 
et al. (2001) suggested that Q « 0.5 is a good approximation 
when the angle dispersion is small (66 < 25°) from numeri- 
cal simulations. From the observations of the molecular con- 
densation in the CO and 13 CO J = 1 — > lines, Takaba et al. 
(1986) estimated an H 2 column density of 2.5 x 10 22 cm~ 2 
and a size of 2.0 pc. The density p can be derived by assuming 
that the line-of-sight size of the dense condensation is similar 
to the lateral size. The FWHM line width of the C ls O J = 1 
— ► line, 2.7 km s^ 1 (Takaba et al. 1986), can be used to es- 
timate Sv\ os . Then the derived field strength is B p w 130 pG. 
The uncertainty in this estimate may be rather large because the 
observed HH 1-2 field is only a part of the L1641 cloud, and 
it should be taken as an order-of-magnitude estimate. The esti- 
mated magnetic field strength of the HH 1-2 region is similar 
to that of other molecular clouds (20-200 p,G) derived using the 
Chandrasekhar-Fermi method (e.g., Andersson & Potter 2005; 
Poidevin & Bastien 2006; Alves et al. 2008). 

An interesting issue is how good near-IR polarimetry is in 
tracing the magnetic field structure of dense clouds. Goodman 
et al. (1995) suggested that the polarizing power of dust grains 
may drop in the dense interior of some dark clouds and that 
near-IR polarization maps of background sources may be un- 
reliable. However, the relevant physics is surprisingly complex 
(Lazarian 2007), and there are observational evidence and the- 
oretical explanations for aligned grains in dense cloud cores 
(Ward-Thompson et al. 2000; Cho & Lazarian 2005). In the 
case of our study of the HH 1-2 field, the polarization degree 
Ph does not show a clear sign of saturation up to H — K s w 1 
(Fig. 7), and the observed polarizations of detected sources may 
well represent the magnetic fields in the HH 1-2 region. 

4.5. P re-Main-Sequence Stars with Infrared Excess 

Sources in group B are PMS stars with infrared excess. In the 
study of global magnetic fields, we excluded group B sources 
because they may have intrinsic polarizations. To verify this 
precaution, their polarization properties may be compared with 
those of group Ar. 

Figure 12 shows the histogram of the polarization position 
angles of group B sources from this work and three bright 
PMS stars from the literature. The distribution is relatively 
widespread, and the peak angle is ill-defined. Probably it is not 
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Fig. 12. — Histogram of near-IR polarization position angles for 
the group B sources and the bright PMS stars (Strom et al. 1985; 
Casali 1995). Note that a 20° bin is used because the sample size 
is small. 



a single-peaked distribution (cf. Tamura & Sato 1989), but the 
sample size is too small to discuss in detail. Quantifying the 
dispersion is difficult, but it is clearly much larger than that of 
group Ar. This distribution implies that these PMS stars have 
significant intrinsic polarizations and that the orientation of the 
intrinsic polarization has no strong correlation with that of the 
global magnetic fields. 

4.6. Outflow Orientations 

The alignment between the magnetic fields of a star-forming 
cloud and the orientation of resulting stars is an interesting 
topic, because it can provide information on the role of mag- 
netic fields during the protostellar collapse and subsequent evo- 
lution of the system. In an ideal scenario of quasi-static isolated 
single star formation, one would expect that the symmetry axis 
of the newly formed (proto)stellar system (such as the rotation 
axis, bipolar outflow axis) would be parallel to the magnetic 
field that threaded the initial dense cloud core (e.g., Galli & Shu 
1993), which would imply that the distribution of outflow ori- 
entations would be similar to that of the polarization angle of 
background sources (with a correction for the averaging effects, 
see below). Strom et al. (1986) found that about 70% of the op- 
tical flows in their imaging survey have the outflow directions 
within 30° of the magnetic field direction of the cloud they be- 
long to. However, in a recent survey of classical T Tauri stars, 
Menard & Duchene (2004) found that the jets/disks around T 
Tauri stars are oriented randomly with respect to the cloud mag- 
netic fields. 

In comparing the distributions of the polarization angle and 
the outflow direction, the dispersion of the polarization angle 
should be corrected by a certain factor, because the averaging 
effects along the line of sight would decrease the dispersion. 
While it is impossible to find an exact correction factor for each 
cloud, numerical simulations suggests that it is in the range of 
0.46-0.51 when the dispersion is small (Ostriker et al. 2001). 
Andersson & Potter (2005) tackled this problem using Monte 
Carlo simulations to analyze the polarimetry of the Southern 
Coalsack cloud and found that a reasonable correction factor 
would be in the range 0.3-1.0. This factor tends to be large 
(close to 1) when the cloud is nearly homogeneous and can 
be small when there are many distinct regions along the line 
of sight. The HH 1-2 region is not completely homogeneous, 
as the molecular cloud has two velocity components along the 
line of sight (see Section 4.3), and also not as complicated as 
the Southern Coalsack cloud. Therefore, a correction factor of 
~0.5 would be a reasonable value. 

For the HH 1-2 region, previously much attention was paid 
to the relation between the direction of the HH 1-2 outflow 
(position angle « 148°) and the orientation of global magnetic 
fields (Strom et al. 1985; WS). Our polarimetry shows that the 
difference between them (^40°) is much larger than the 11° 
dispersion of the polarization angles. Therefore, there is a cer- 
tain degree of misalignment. To make a more meaningful com- 
parison, a list of optical jets in the HH 1-2 field with known 
flow directions was compiled (Table 3). The driving sources of 
these outflows are YSOs. Figure 13 shows the distribution of 
the outflow orientations. Considering the distribution of the po- 
larization angle and the correction factor mentioned above, the 
expected distribution of the outflow direction would show a sin- 
gle peak at ~1 10° with a dispersion of ^20°. Surprisingly, the 
position angle of the outflow shown in Figure 13 seems to be 
nearly random, i.e., the distribution is almost uniform, which 
makes a stark contrast to the distribution of polarization angles. 
Therefore, we suggest that many protostars may become disori- 
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Fig. 13. — Histogram of outflow orientations of optical jets in 
the HH 1-2 region. 



ented during the star forming process. 

The random orientation of YSOs suggests that the global 
magnetic field may lose control of a collapsing protostellar core 
at a certain small scale or at a certain stage of evolution. Several 
mechanisms may be working to cause the disorientation. First, 
the dynamical interaction during a binary (or multiple star) for- 
mation process can produce a misaligned system (e.g., Bonnell 
et al. 1992). In fact, misaligned binaries in the Class II phase 
are not unusual (Monin et al. 2007). Among the outflow driv- 
ing sources in the HH 1-2 region, V380 Ori and HH 1-2 VLA 
1/2 are examples to the point. V380 Ori is a 0'.'15 (-60 AU) bi- 
nary and drives HH 35 and HH 148, and these two outflows are 
almost perpendicular to each other (Strom et al. 1986; Leinert 
et al. 1997). HH 1-2 VLA 1/2 is a 3" (-1200 AU) binary and 
drives HH 1-2 and HH 144-145, and the projected angle be- 
tween these flows is —70° (Reipurth et al. 1993). Second, even 
in a single-star formation, an outflow may not be aligned with 
the magnetic field of the surrounding cloud if the magnetic field 
is too weak (Matsumoto et al. 2006). There are other possible 
mechanisms (Menard & Duchene 2004). 

Considering that there are Class sources showing misalign- 
ments, the disorientation may happen very early. For example, 
HH 1-2 VLA 1 is a Class source (Chini et al. 1997). Exam- 
ples in other star-forming regions include HH 24 MMS, which 
seems to have two accretion disks with a ~45° difference in 
projected orientation (Kang et al. 2008), and NGC 1333 IRAS 
2, which seems to drive two outflows almost perpendicular to 
each other (Hodapp & Ladd 1995). 

4. 7. Wavelength Dependence of Interstellar Polarization 

Though the selected point-like sources have no detectable 
nebulosity, we cannot rule out the possibility of the presence of 
unresolved reflection nebulae in some of the sources. Indeed, 
Casali (1995) showed that the polarization by scattering is im- 
portant in some of the sources in LI 641, especially the sources 
with low extinction. To discriminate between the contributions 
from the dichroic extinction and from the scattering, the wave- 
length dependence of polarization can be measured by calculat- 
ing the ratio of polarization degrees. In the infrared wavelength 
range, the polarization by dichroic extinction decreases with 
wavelength while the polarization by scattering is not a strong 
function of wavelength (Whittet et al. 1992; Casali 1995). 

Figure 14 shows the Pj/Pk s and Ph/Pk s ratios for the 
sources in the HH 1-2 field. It is very clear that groups AO and 
Ar show very different behavior. Within group Ar, the ratios 
are reasonably constant, and the weighted average values are 
Pj/Pk s =2.18 and Ph/Pk s = 1-51 with standard deviations 
of 0.10 and 0.07, respectively. These values are consistent with 
those for the whole L1641 area (2.5 and 1.4 with 2a scatters of 
0.7 and 0.2, respectively; Casali 1995), which is also consistent 
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Fig. 14. — The ratio of J to K s polarizations {top panel) and 
the ratio of H to K s polarizations (bottom panel) against the 
H — K s color for the 26 sources detected in the J, H, and K s 
bands. Triangles: sources in group AO. Source 30 has unusually 
low ratios. Filled circles: sources in group Ar. Squares: sources 
in group B. Dotted lines: weighted average for group Ar. 

with the empirical relation P oc with (i = 1.6-2.0 (Whittet 
1992). Therefore, the polarization of group Ar sources can be 
very well explained by dichroic extinction. 

In contrast, the Pj/Pk 3 and Ph/Pk 3 ratios of group AO 
sources are near unity, which is significantly different from 
those of group Ar sources. Therefore, the polarization mech- 
anism for low extinction sources may be dominated by the 
circumstellar scattering (see more discussions in Section 5.4 
of Casali 1995). Source 30 shows unusually low ratios (Fig. 
14) because its K s band polarization degree is high while the 
Pj I Ph ratio is near one as expected. Its K s band polarization 
angle is also quite different from those of J and H bands. It is 
not clear what caused this peculiar behavior. 

5. SUMMARY 

We conducted a deep and wide-field J-H-K s imaging po- 
larimetry toward a 8' x 8' region around HH 1-2 in the star- 
forming cloud LI 641. The main results in this study are sum- 
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marized as follows. 
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Table 1 

Photometry of Point-like Sources in the HH 1-2 Field 
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Position 
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15.78 


Ar 


43 


5 36 09.81 


-6 44 09.3 


17.39 


15.86 


15.14 


Ar 


44 


5 36 07.19 


-6 44 08.8 


19.04 


16.90 


15.91 


Ar 


45 


5 36 36.38 


-6 44 08.1 


18.99 


17.84 


16.95 


B 


46 


5 36 13.28 


_6 44 02.4 


18.08 


16.14 


15.18 


Ar 


47 


5 36 11.37 


-6 44 00.1 


16.69 


15.23 


14.55 


Ar 


48 


5 36 13.45 


-6 43 54.7 


18.00 


16.88 


16.22 


B 


49 


5 36 29.06 


-6 43 51.5 


17.83 


15.98 


15.11 


Ar 


50 


5 36 26.77 


-6 43 43.4 


13.84 


11.94 


11.07 


Ar 


51 


5 36 12.69 


-6 43 34.0 


15.81 


14.35 


13.64 


Ar 
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TABLE 1 — Continued 



Source 


Position 


J 


H 


K s 


Group a 




Ct T9nnn n 


o T9nnn n 


(mag) 


(mag) 


(mag) 




52 


5 36 15.27 


-6 43 30.8 


19.22 


17.02 


1 c nt 

16.01 


A „ 

Ar 


53 


5 36 36.63 


-6 43 23.2 


1 C 11 

15.72 


1/1 A^l 

14.42. 


1 1 oo 

13.88 


A „ 

Ar 


C A 

54 


5 36 34.28 


-6 43 23.2 


13.1/ 


12.51 


12.33 


AO 


55 


5 36 27.59 


-6 43 22.2 


18.54 


16.69 


15.92 


Ar 


56 


5 36 30.18 


-6 43 20.4 


17.19 


1 C HH 

15.77 


1 C 1 A 

15.14 


Ar 


57 


5 36 20 52 


-6 43 18.1 


16.91 


16.39 


16.10 


AO 


58 


5 36 07.34 


-6 43 07.6 


12.89 


11.88 


11.42 


Ar 


59 


5 36 35.42 


-6 43 07.5 


18.69 


17.16 


16.50 


Ar 


60 


5 36 25.94 


-6 43 02.2 


13.85 


13.01 


12.55 


Ar 


61 


5 36 35.76 


-6 42 49.9 


13.48 


12.38 


11.93 


Ar 


62 


5 36 36.87 


_6 42 49.4 


19.10 


17.02 


16.07 


Ar 


63 


5 36 07.61 


-6 42 46.6 


18.14 


16.48 


15.73 


Ar 


64 


5 36 30.23 


-6 42 46.1 


13.31 


11.81 


11.09 


Ar 


65 


5 36 12.83 


-6 42 34.6 


19.06 


17.26 


16.30 


Ar 


66 


5 36 35.14 


-6 42 18.6 


18.25 


16.66 


16.00 


Ar 


67 


5 36 30.53 


-6 42 03.1 


14.99 


14.44 


14.16 


AO 


68 


5 36 20.82 


-6 41 56.3 


19.52 


17.53 


16.71 


Ar 


69 


5 36 32.13 


-6 41 51.6 


16.78 


15.54 


15.03 


Ar 


70 


5 36 17.49 


-6 41 46.1 


18.28 


17.08 


16.42 


Ar 


71 


5 36 21.96 


-6 41 42.0 


12.81 


12.22 


11.93 


AO 


72 


5 36 35.81 


-6 41 41.3 


16.95 


16.38 


16.05 


AO 


73 


5 36 18.80 


-6 41 28.9 


18.06 


16.16 


15.30 


Ar 


74 


5 36 35.37 


-6 41 29.2 


16.08 


14.84 


14.33 


Ar 


75 


5 36 19.23 


-6 41 18.2 


16.69 


15.16 


14.46 


Ar 


76 


5 36 31.52 


-6 41 13.6 


12.72 


11.98 


11.78 


AO 



NOTES. — Units of right ascension are hours, minutes, and seconds, and units 
of declination are degrees, arcminutes, and arcseconds. Positions are from the 
J-H-K s image (Fig. 1). 

a Classification based on a color-color diagram (see Section 4.2). 
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Table 2 

polarimetry of point-like sources in the hh 1-2 field 



— . 

Source 




p.. 

*H 


p„ 

*Ks 


ft T 


ft „ 


ft „ 






\ /(! ) 


V >°) 


v t 


I ) 


\ ) 


i 


< 12.2 


< 6.2 


< 20.4 








2 


< 9.6 


3.6 ± 1.0 


< 9.5 




113.6 ± 8.0 




3 


40.5 ± 8.5 


< 10.8 


< 32.8 


109.2 ± 5.9 






4 


2.74 ± 0.08 


2.27 ± 0.04 


1.55 ± 0.09 


97.7 ± 0.8 


95.7 ± 0.5 


83.3 ± 1.6 


5 


< 9.4 


< 4.4 


< 14.4 








6 


31.3 ± 10.2 


5.7 ± 1.8 


< 12.7 


51.1 ± 8.9 


111.8 ± 8.5 




7 


2.84 ± 0.17 


2.23 ± 0.09 


1.1 ±0.3 


121.3 ± 1.7 


121.4 ± 1.1 


110.8 ± 7.0 


8 


0.30 ± 0.04 


0.44 ± 0.03 


0.42 ± 0.08 


108.7 ± 3.7 


120.8 ± 2.0 


94.9 ± 5.4 


9 


0.33 ± 0.07 


0.70 ± 0.05 


< 0.5 


116.2 ±6.3 


109.7 ± 2.1 




10 


10.97 ± 0.06 


9.40 ± 0.03 


7.30 ± 0.04 


148.6 ±0.1 


147.7 ± 0.1 


147.9 ±0.1 


11 


4.3 ± 0.7 


3.4 ±0.3 


< 4.3 


143.6 ± 4.8 


127.9 ± 2.9 




12 


2.4 ± 0.6 


2.2 ±0.3 


< 2.7 


119.8 ± 7.1 


116.0 ± 3.6 




13 


< 13.6 


< 5.0 


< 13.2 








14 


6.2 ± 0.9 


4.1 ±0.5 


< 5.2 


118.6 ±4.0 


116.9 ± 3.6 




15 


1.8 ± 0.6 


2.1 ±0.2 


< 2.1 


122.2 ± 8.8 


125.3 ± 3.2 




16 


< 8.4 


2.9 ±0.8 


< 5.9 




121.7 ± 7.7 




17 


7.2 ± 2.1 


3.1 ±0.3 


2.8 ±0.6 


127.8 ± 7.9 


133.5 ± 3.2 


119.8 ± 5.6 


18 


4.5 ± 0.5 


3.3 ±0.2 


3.2 ±0.7 


115.5 ± 3.0 


118.7 ±2.1 


123.8 ± 6.3 


19 


< 30.2 


< 7.0 


< 19.2 








20 


< 10.2 


7.3 ± 1.2 


< 8.8 




132.1 ±4.7 




21 


7.7 ± 0.6 


6.07 ± 0.15 


3.9 ±0.2 


20.3 ± 2.2 


21.0 ±0.7 


14.0 ± 1.7 


22 


< 2.7 


2.8 ±0.5 


< 6.4 




118.2 ±5.3 




23 


< 1.2 


< 0.9 


< 3.4 








24 


< 8.3 


4.4 ±0.4 


5.0 ±0.8 




84.0 ± 2.6 


82.6 ± 4.5 


25 


6.9 ± 0.7 


5.00 ±0.10 


3.21 ± 0.16 


100.3 ± 2.9 


106.6 ± 0.6 


108.5 ± 1.4 


26 


0.67 ± 0.17 


0.96 ±0.12 


< 1.1 


110.0 ±7.1 


126.0 ± 3.6 




27 


5.89 ± 0.18 


4.03 ± 0.03 


2.59 ± 0.03 


114.6 ±0.9 


113.5 ±0.2 


112.5 ±0.3 


28 


7.3 ± 1.6 


4.5 ±0.5 


< 3.5 


118.9 ±5.9 


127.4 ± 2.9 




29 


0.51 ± 0.05 


0.74 ± 0.04 


0.53 ± 0.09 


118.9 ±3.0 


112.8 ± 1.5 


99.6 ± 5.0 


30 


0.62 ± 0.20 


0.65 ± 0.16 


2.4 ±0.6 


107.3 ± 8.7 


123.3 ± 6.7 


65.8 ± 6.6 


31 


4.8 ± 1.3 


< 2.0 


7.8 ±2.3 


115.5 ± 7.5 




137.9 ±8.1 


32 


13.6 ± 2.5 


5.2 ±0.4 


3.5 ±0.8 


79.0 ±5.1 


82.6 ± 2.2 


72.3 ± 6.3 


33 


< 9.6 


5.9 ±0.9 


< 5.6 




114.2 ±4.4 




34 


6.5 ± 0.5 


5.2 ±0.2 


4.1 ±0.5 


126.0 ± 2.2 


115.1 ± 1.1 


130.4 ±3.3 


35 


0.71 ± 0.05 


0.83 ± 0.03 


0.65 ± 0.06 


71.6 ±2.1 


71.3 ± 1.1 


70.9 ± 2.8 


36 


1.49 ± 0.10 


1.12 ± 0.04 


0.68 ± 0.08 


118.0 ±2.0 


109.5 ± 1.1 


109.1 ± 3.3 


37 


< 3.5 


< 1.0 


1.7 ±0.5 






118.8 ± 8.1 


38 


2.9 ± 0.8 


2.3 ±0.2 


1.8 ±0.4 


124.5 ± 8.1 


111.6 ±2.7 


109.0 ± 6.3 


39 


3.9 ± 1.0 


3.2 ±0.5 


< 5.7 


107.8 ± 6.7 


113.1 ±4.5 




40 


1.59 ± 0.05 


1.38 ±0.03 


0.93 ± 0.03 


132.8 ±0.9 


128.0 ± 0.5 


126.7 ± 0.9 


41 


0.57 ± 0.05 


0.65 ± 0.04 


0.48 ± 0.09 


115.9 ±2.7 


114.8 ± 1.6 


112.8 ± 5.4 


42 


< 8.2 


5.0 ± 1.2 


< 7.5 




114.7 ±6.4 




43 


< 4.3 


2.4 ±0.5 


< 3.8 




107.8 ± 5.3 




44 


< 26.1 


5.6 ± 1.7 


< 9.5 




126.2 ± 8.2 




45 


< 15.0 


< 10.0 


< 20.7 








46 


10.9 ± 2.5 


5.4 ±0.7 


< 4.1 


82.9 ± 6.3 


107.1 ± 3.9 




47 


3.5 ± 0.7 


2.0 ±0.3 


< 2.3 


84.2 ± 5.9 


113.6 ± 3.7 




48 


17.8 ± 2.5 


8.8 ± 1.3 


< 10.1 


149.5 ± 4.0 


133.9 ±4.1 




49 


< 7.5 


5.2 ±0.6 


< 3.8 




104.6 ± 3.2 




50 


1.64 ± 0.11 


1.09 ±0.04 


0.56 ± 0.04 


90.3 ± 1.8 


91.8 ± 1.0 


109.0 ± 2.0 


51 


4.7 ± 0.4 


3.03 ±0.14 


1.8 ±0.3 


86.9 ±2.1 


95.3 ± 1.3 


99.2 ± 4.9 


52 


< 21.9 


< 4.2 


< 8.7 








53 


3.6 ± 0.3 


2.56 ±0.16 


1.7 ±0.4 


131.5 ± 2.7 


122.8 ± 1.8 


137.9 ±6.5 


54 


0.66 ± 0.05 


0.70 ± 0.04 


0.63 ± 0.10 


125.6 ± 2.4 


120.7 ± 1.4 


110.7 ±4.6 


55 


< 17.8 


8.3 ± 1.4 


10.3 ±2.8 




106.4 ± 4.8 


128.4 ±7.6 


56 


< 3.7 


3.4 ±0.5 


< 3.7 




110.1 ±4.1 




57 


< 3.1 


2.6 ±0.8 


< 9.4 




104.5 ± 8.1 




58 


4.14 ± 0.05 


2.54 ± 0.03 


1.39 ± 0.06 


92.9 ± 0.3 


95.2 ± 0.3 


103.5 ± 1.2 


59 


< 13.7 


< 4.9 


< 13.5 








60 


< 3.6 


< 2.8 


< 2.5 








61 


2.31 ± 0.07 


1.70 ± 0.04 


1.08 ± 0.07 


130.5 ± 0.8 


125.4 ± 0.6 


114.5 ± 1.9 


62 


< 17.7 


7.0 ± 1.5 


< 9.8 




135.6 ±5.9 




63 


< 10.6 


< 2.8 


< 6.6 








64 


1.47 ± 0.06 


1.21 ±0.03 


0.81 ± 0.04 


126.7 ± 1.1 


129.3 ± 0.6 


138.5 ± 1.3 


65 


< 17.3 


< 4.7 


< 11.3 








66 


13.7 ± 3.1 


< 3.8 


< 8.8 


120.3 ± 6.3 
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TABLE 2 — Continued 



Source 


Pj 


Ph 


Pk s 


9j 


Oh 






(%) 


(%) 


(%) 


(°) 


(°) 


(°) 


67 


< 0.5 


0.52 ±0.15 


2.0 ±0.5 




98.6 ± 8.1 


124.5 ± 7.3 


68 


< 33.9 


< 7.5 


< 16.2 








69 


7.2 ± 0.8 


4.3 ±0.4 


< 3.4 


115.7 ± 3.2 


118.6 ±2.6 




70 


9.9 ± 3.1 


< 4.6 


13.1 ±4.3 


35.6 ± 8.6 




100.9 ± 8.9 


71 


0.43 ± 0.05 


0.62 ± 0.03 


0.45 ± 0.07 


135.8 ± 3.0 


134.6 ± 1.4 


134.3 ± 4.6 


72 


< 4.0 


< 2.5 


< 9.2 








73 


< 10.0 


5.1 ±0.8 


< 5.2 




101.3 ± 4.2 




74 


2.7 ± 0.7 


3.5 ±0.2 


< 2.2 


120.4 ± 7.5 


120.9 ± 2.0 




75 


8.8 ± 1.4 


4.6 ±0.4 


4.6 ± 1.0 


88.1 ± 4.3 


101.0 ± 2.6 


123.3 ± 6.3 


76 


0.58 ± 0.06 


0.53 ± 0.03 


< 0.34 


110.0 ± 2.8 


105.0 ± 1.8 




Notes. - 


—For sources with P/SP < 3, the 3SP upper limits are listed. 







Table 3 

Optical Jets in the HH 1-2 Field 



Object 


P.A. a 


Driving Source 


References 


HH 1-2 


148° 


VLA 1 


Pravdo et al. 1985 


HH 35 


149° 


V380 Ori 


Strom et al. 1986 


HH 144-145 


82° 


VLA 2 


Reipurth et al. 1993 


HH 146 


6° 


VLA 4 


Reipurth et al. 1993 


HH 147 


50° 


N 3 SK50 


Eisloffel et al. 1994 


HH 148 


56° 


V380 Ori 


Strom et al. 1986 


SMZ61 


171° 


VLA 3 


Stanke et al. 2002 



a Position angle of the outflow axis. Typical uncertainty is 10°. 
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